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We study the absorbing phase transitions in coupled directed percolation (DP) processes with A''- 
species particles in one dimension. The interspecies coupling is linear, bidirectional, and excitatory. 
We find that the presence of a spontaneous annihilation process A — > is essential in stabilizing the 
absorbing phase (vacuum). In the coupled contact processes, the vacuum is stable and the system 
exhibits DP type transitions, regardless of the coupling strength, for all A'^. However, in the coupled 
branching annihilation random walks with one offspring (BAW), where particle annihilations occur 
only through binary diffusion processes A + A — > 0, the vacuum becomes unstable with respect to 
an arbitrarily small branching rate in a sufficiently strong coupling regime for A' > 3. The N — 2 
BAW exhibits the DP type transition for any coupling strength, but the inclusion of interspecies 
hard core (HC) interaction makes the vacuum unstable again and the system is always active in 
a strong coupling regime. Critical behavior near the zero branching point is characterized by the 
scaling exponents, (5 = v± = 1/2 and v\\ = 1, regardless of the presence of HC interaction. We also 
discuss the effects of the asymmetric coupling. 
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I. INTRODUCTION 

Nonequilibrium systems with absorbing (trapped) 
states have been extensively studied in recent years, be- 
cause of their wide applicability to various phenomena in 
natural science as well as in social and economical science 
P, 13 ■ The absorbing transition between an active phase 
into an absorbing phase can be regarded as one of the 
simplest and natural extensions of the well-established 
equilibrium phase transition to nonequilibrium systems. 

Critical behaviors near an absorbing transition are cat- 
egorized into a few universality classes characterized by 
the symmetry between absorbing states and/or the con- 
servation in dynamics |^ ^ |^ . The most prominent uni- 
versality class is the directed percolation (DP) class 
which involves typically a single absorbing state without 
any conservation in dynamics. Most of absorbing-type 
nonequilibrium models belong to the DP class. Only a 
small number of models form different classes such as the 
directed Ising (DI) class involving two equivalent classes 
of absorbing states 0, 0, 13 and the parity-conserving 
(PC) class with mod (2) conservation of the total particle 
number . These two classes coincide in one dimension 
by identifying a domain wall in the Dl-type models as 
a particle in the PC-type models. In higher dimensions, 
both types of models are always active (no absorbing 
phase) except at a trivial point (annihilation fixed point) 
and their critical behavior is described by the mean field 
theory. None of the models studied so far with higher 
symmetries than the Ising symmetry (for example. Potts 
symmetry) stabilizes an absorbing phase even in one di- 
mension. 
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Recently, a variety of coupled systems have been in- 
vestigated extensively. Janssen ^ studied coupled DP 
processes with bilinear and bidirectional interspecies cou- 
plings in the framework of bosonic field theory, where 
no other critical phenomena were found than the DP. 
Tauber et al [lOj studied linearly and unidirectionally 
coupled DP processes, where a series of new multicrit- 
ical phenomena was observed. Coupled PC (DI) pro- 
cesses have been also studied where the absorb- 
ing phase become unstable with respect to an arbitrarily 
small branching rate even in one dimension. More inter- 
estingly, the critical behavior near the annihilation fixed 
point depends on details of particle dynamics such as the 
presence of interspecies hard core (HC) interaction and 
the branching method [T^ IT^j l . 

Stochastic models for linearly and bidirectionally cou- 
pled DP processes have been also studied through two- 
species branching annihilating random walks with one 
offspring (BAWi) in one dimension It was found 

that the HC interaction is crucial in a strong coupling 
regime, where the absorbing phase (vacuum) becomes 
unstable and the system is always active except at the 
annihilation fixed point of zero branching rate. Critical 
behavior near the annihilation fixed point was conjec- 
tured via an analytic argument and confirmed by numer- 
ical simulations 14]. 

In this paper, we studied the iV-species BAWi models 
and also the A^-species contact processes (CP) with and 
without interspecies HC interaction. The single-species 
BAWi and CP both belong to the DP universality class, 
but their upper critical dimensions are different (due = 2 
and 4) HEl El- The essential difference is the ab- 
sence of spontaneous annihilation process A ^ in the 
BAWi, where particle annihilations occur only in pairs 
through binary diffusion processes A + A ^ This an- 
nihilation random walk nature of the BAWi makes the 
vacuum unstable for d > 2. In the coupled systems, one 
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can expect a similar scenario that the stabiHty of the vac- 
uum also depends crucially on the presence of the sponta- 
neous annihilation process. Furthermore, as increases, 
the vacuum becomes more unstable in the 7V-B AWi , be- 
cause the annihilation process of only the same species 
of particles is allowed. A sufficiently strong interspecies 
coupling may wash away completely the absorbing phase 
for large enough N even in one dimension, even without 
any HC interaction. In this paper, we address this vac- 
uum stability question in the strong coupling regime and 
measure the threshold value of N with and without HC 
interaction. 

Outline of this paper is as follows. In Sec. II, the N- 
species CP models are introduced. We find that, for any 
N, the vacuum is stable for a low branching rate and 
the DP-type phase transition is observed regardless of 
the presence of HC interaction. In Sec. HI, we study 
the A^-species BAWi models by numerical simulations. 
Wc find that, for iV > 3, the vacuum becomes unstable 
with respect to an arbitrarily small branching rate in a 
sufficiently strong coupling regime. The critical behav- 
ior near the zero branching point is explored in terms of 
the scaling exponents. In a weak coupling regime, the 
conventional DP transition into a stable vacuum is ob- 
served. The interspecies HC interaction only shifts the 
DP critical points and the threshold value of the inter- 
species coupling for complete disappearance of a stable 
vacuum. The N = 2 case is special. Without HC in- 
teraction, one can always find the DP transition into a 
stable vacuum even at the full coupling strength, like in 
the iV-species CP models. However, the inclusion of HC 
interaction makes the vacuum unstable again in a strong 
coupling regime. We also discuss the effect of the direc- 
tional asymmetry in the bidirectional couplings on the 
phase diagram. Finally, we conclude in Sec. IV. 

II. COUPLED CONTACT PROCESSES 

The A^-species coupled contact process (iV-CP) is de- 
fined by the following evolution rules: (1) Each parti- 
cle annihilates spontaneously with probability p or (2) 
creates a particle of the same species with probability 
cr(l — p) or a particle of the different species with prob- 
ability (t'{1 ~ p)/{N — 1) in its neighborhood (branching 
process) as 

Ai — > withp, 

A, — > A, + Ai with cr(l -p), (1) 
A, — > A, + Aj with cr'(l -p)/(A^- 1), 

where a' ^ 1 — a, i ^ j, and i — I, . . . , N . Any branch- 
ing attempt is rejected if it would result in a multiple 
occupation by the same species particles at a site. With 
interspecies HC interaction, a multiple occupation by dif- 
ferent species particles is also forbidden. Each species is 
symmetrically coupled with other species by the inter- 
species coupling strength a' . For a' — 0, all species are 
completely decoupled. 
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FIG. 1: The a -p phase diagram for the A^-CF, (a) A' = 2 
and (b) TV = 3. Filled and open circles correspond to critical 
points with and without HC interactions, respectively. Lines 
between data points are guides to the eyes only. 

It is trivial to show that the A^-CP models with HC 
interaction become identical to the single-species CP 
model. Each site can be occupied by only one particle, 
regardless of its species. Any configuration and any dy- 
namic process in the N-CP can be exactly mapped on its 
corresponding configuration and dynamic process of the 
single-species CP by simply ignoring the particle species. 

Without HC interaction, a multiple occupation by dif- 
ferent species particles is allowed. One can easily expect 
that the active phase expands until the multiple occupa- 
tion is maximized. For N = 2, there exists a simple dual- 
ity between two species {a ^ a') and the maximum point 
is located at cr = ct' = 1/2, see Fig.[T](a). For A^ > 3, the 
maximum point is expected to be at ct = 1/N, where the 
branching symmetry is perfect and the maximum mix- 
ing is expected, see Fig. ^(b). From this reasoning, we 
expect that complete disappearance of the vacuum does 
not occur at any coupling strength a for any finite A^. Of 
course, all absorbing critical phenomena should belong 
to the DP class. 



III. COUPLED BAWi 

The A^-species coupled branching annihilating random 
walks with one offspring (A^-BAWi), is defined by the 
following evolution rules: (1) Each particle hops to a 
nearest-neighboring site with probability p or (2) creates 
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a particle of the same species with probabihty (t(1 — p) or 
a particle of the different species with probability cr'(l — 
p)/{N — 1) in its neighborhood (branching process) as 

AS < — > 9Ai with p, 

Ai — > Ai + A, withcr(l-p), (2) 
Ai — > Ai + Aj with cr'(l -p)/(iV - 1), 

where cr' 1 — cr and i ^ j. If two identical particles hap- 
pen to be on the same site, both particles immediately 
annihilate each other {Ai + Ai — > 0). With interspecies 
HC interaction, any hopping or branching attempt is re- 
jected if it would result in a multiple occupation of par- 
ticles, regardless of their species. The above model has 
a permutational symmetry between species. Later, we 
will also study an asymmetrically coupled system which 
breaks the permutational symmetry. 

A. N = 2 

The two-species BAWi model has been investigated 
previously for general a In this paper, we present 

the numerically improved phase diagram in Fig. [5] (a) 
and briefly summarize the results of |l4j for comparison 
to those for > 3 and also for self-containedness. With- 
out HC interaction, the system always exhibits the DP- 
type absorbing transition for all cr, including the case 
of the maximum coupling strength {a' — 1 or cr = 0). 
However, with HC interaction, there exists a strong cou- 
pling regime (ct' > a'*) where the vacuum becomes un- 
stable with respect to an arbitrarily small branching rate 
[p = 1~) and completely disappears in the phase dia- 
gram except at the annihilation fixed point {p — 1). The 
threshold value for the strong coupling regime is numer- 
ically estimated as cr'* ~ 0.95. In the weak coupling 
regime (cr' < cr'*), we observe the conventional DP tran- 
sition into the vacuum. 

It is not surprising to see that the active phase ex- 
pands as ct' increases, because a collision chance of the 
same species particles decreases. However, the effect of 
HC interaction is rather tricky. For small cr', the system 
tends to form large domains of the same species parti- 
cles. The HC interaction induces an effective diffusion 
bias directed to the domain center, which accelerates the 
pair annihilation process. Therefore the system becomes 
less active with HC interaction. For large cr', the situ- 
ation is reversed. The system prefers locally heteroge- 
neous configurations and the HC interaction reduces a 
chance of binary collision of the same species particles. 
Thus, in this case, the system becomes more active with 
HC interaction. This explains why two critical lines with 
and without HC interaction should cross each other as in 
Fig.H 

Consider the ct = line, where a particle cannot cre- 
ate a particle of the same species. As a result, a single 
particle cannot be annihilated by a single branching and 
diffusion process: A ^ A + A ^ $. It needs at least three 
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FIG. 2: The a-p phase diagram of the iV-BAWi , (a.) N = 2 
and (b) A'' = 3. Filled and open circles correspond to critical 
points with and without HC interactions, respectively. Lines 
between data points are guides to the eyes only. 

branching processes such as A ^ AB ABA ABAB 
or ABBA. This four particle state can turn into vacuum 
via successive pair annihilations only if diffusion across a 
different species particle is allowed (no HC interaction). 
Hence, the vacuum can be stable in a low branching (high 
diffusion) regime. However, in the presence of HC in- 
teraction, the ordered AB pairs can not be annihilated 
by diffusions. Therefore, a single particle has a nonzero 
probability to survive asymptotically and the vacuum is 
unstable with respect to an arbitrarily small branching 
rate. 

The critical behavior near the annihilation fixed point 
in the strong coupling regime is characterized by a set of 
the scaling exponents 

/3= 1/2, = 1/2, = 1, (3) 

where the exponents (3, and characterize the scal- 
ing behavior of the steady-state particle density ps, the 
correlation length ^, and the relaxation time r, respec- 
tively. 

The exponent ratios, P/h'^ = 1 and z = v\\lv^ = 2, 
originate from the ordinary diffusion nature at the an- 
nihilation fixed point in one dimension 12] . There is 
only one independent exponent /3, of which the value 
can be extracted by a simple argument 0|. Consider 
a particle A created by a particle B. This branching 
process increases the particle density with the time scale 
Tb (1 —p)~^. Near the annihilation fixed point of zero 



4 



branching rate, the offspring A would be annihilated by 
colliding with an independent A via diffusion. The time 
scale for this process is governed by ordinary diffusion: 
Td ~ £^ where the mean distance between particles £ is 
order of the inverse of the particle density p~^. Balanc- 
ing these two time scales, we can expect the steady-state 
particle density to scale as ps ~ (1 — p)^ with f3 = 1/2. 

This argument is quite general, so it should apply to 
many other models exhibiting a critical behavior near 
the annihilation fixed point. Moreover, the HC interac- 
tion does not matter in this argument. So, one can expect 
that the A^-species BAWi should belong to the same class 
(/3 = 1/2) for any N >2, regardless of the presence of HC 
interaction, which will be confirmed numerially in next 
subsection. Moreover, the A^-species BAW2 also belongs 
to the same class when the HC interaction is present and 
the branching method is static: Two offspring are divided 
by their parent, 0B0 ABA. It is clear that our argu- 
ment applies to this model, due to the HC interaction. 
When the branching method is dynamic or there is no 
HC interaction, A^-BAW2 belongs to different universal- 
ity classes [T^. 



B. N>3 

First, consider the = 3 case without HC interac- 
tion. To map out the phase diagram, we perform the 
defect- type dynamic Monte Carlo simulations, starting 
with a single particle. We measure the survival proba- 
bility P{t), the number of particles N{t), the mean dis- 
tance of spreading R{t). At criticality, these quantities 
scale algebraically in the l ong time limit as P{t) ^ t^^ , 



N{t) ~ f and R{t) ~ i^/^ 

By inspecting the curvature of effective exponents de- 
fined as 



6{t)^\og[P{t)/P{t/m)]/\ogm, 



(4) 



with arbitrary m, and similarly for ri(t) and l/z{t), we 
estimate the values of the critical hopping probability Pc 
and the dynamic exponents (5, rj and z for various values 
of a. The a — p phase diagram is shown in Fig. [21 (b). 
For small a, the absorbing phase (vacuum) completely 
disappear. This result is rather unexpected because, for 
N = 2, the vacuum is always stable for a low branch- 
ing rate without HC interaction. Before going into de- 
tailed discussion on this vacuum instability without HC 
interaction, we present numerical results. The threshold 
value of the vacuum instability is estimated as a* ~ 0.125 
(a'* ~ 0.875). For a > a*, we observe that the sys- 
tem undergoes an absorbing transition into vacuum. As 
expected, we find the DP critical exponents along both 
paths of constant p and constant cr lines. 

To identify the scaling behavior near the annihilation 
fixed point {p = 1) for a < a*, we analyze the finite-size 
effects on the steady-state particle density ps ■ Using the 
finite-size scaling theory on ps Tsj 
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FIG. 3: Data collapse of psL''^''"^ against AL^^"^ with 
l3/v_L ~ 1 and v±_ — 1/2 for system size L = 2^, ■ ■ ■ , 2^ for 
= 3 (a) without HC interaction at cr = (b) with HC 
interaction at o — 0.2. 



with A = Pc — p and system size L, the value of iy± is 
determined by collapsing the data of ps with (3/i'± = 1. 
We measure ps in the steady state, averaged over 5 x 
lO'^ ^ 5 X 10^ samples for several values ofA(5xlO^*'^ 
0.05) and L (2^ - 2^). 

We estimate = 0.50(4) at cr = 0. Fig|3| (a) shows 
very good data collapse of ps with = 0.50. These 
numerical results confirm our previous argument for the 
universality class near the annihilation fixed point. 

To understand the vacuum instability for small a for 
A^ = 3 without HC interaction, we again consider the 
cr = line, where a particle can create a particle of dif- 
ferent species only. As explained before, there needs a 
sequence of branching processes for a single particle anni- 
hilation. For N = 2 without HC interaction, the process 
of A^ AB ABA ABAB ^ ^ is possible. 
In this case, AA or BB pairs will quickly go away at a 
very high diffusion rate, so one may consider only AB 
pairs, which cannot be annihilated by itself. The last 
two reaction events can be regarded as a collision and 
annihilation event of two AB pairs and the previous re- 
action events as a creation event of a AB pair out of a 
parent AB pair. This process is exactly equivalent to the 
A^ = 1 BAWi model by identifying a AB pair as a par- 
ticle. Therefore the N = 2 case without HC interaction 
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can have a stable vacuum like in the A'^ = 1 model. The 
iV = 3 case is different. Three different pairs {AB, BC, 
CA) are possible. Any pair can branch any other pair, 
e.g. AB ACBA ^ CB ^ CBBA CA. So, the 
N — 3 case can not reduce to the = 1 model, in con- 
trast to the N — 2 case. It seems that this multispecies 
character is the crucial element for vacuum instability in 
BAW-type binary diffusion-annihilation models. 

Interspecies HC interaction should destabilize the vac- 
uum more easily. We expect that the threshold of vac- 
uum instability should be lowered with HC interaction. 
In Fig. 121 (b), the phase diagram for iV = 3 with HC 
interaction is presented. We estimate a* ~ 0.25. To es- 
timate i>± for a < a* , we also try the data collapse of Ps ■ 
In Fig. 13 (b), we present the data at a = 0.2 with the 
exponent values of v± = 0.50 and I3/va_ = 1. Our esti- 
mation is v±^ = 0.50(2). Again, we find that our simple 
argument also applies to this case. 

We also perform defect-type Monte Carlo simulations 
for the = 4 case. We estimate cr* = 0.22(2) without 
HC interaction and a* — 0.328(2) with HC interaction. 
The exponent value is estimated as v±^ ~ 0.52(2) at u = 
0.15 without HC interaction and v^^ 2± 0.52(2) at u = 0.3 
with HC interaction. 



C. asymmetric coupling 

We study the effect of broken permutational symmetry 
on the phase diagram. For convenience, we consider the 
= 3 case at cr = (ct' = 1) without HC interaction 
only. The branching process in Eq. |(5J) is modified as 
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(6) 



where < q < 1. 

At q — 0, the species C is completely suppressed and 
the model becomes identical to the N — 2 symmetric one 
at cr = 0. At q = 1/2, all three species are equivalent and 
the A^ = 3 symmetric model is recovered. Our results on 
the symmetric A'^-BAWi models in previous subsections 
indicate that the vacuum is stable in high diffusion regime 
at q = 0, but becomes unstable completely at q — 1/2. 
For q < 1/2. the species C is suppressed in comparison 
to the other two species A and B. Therefore, the density 
of the third species C should be proportional to q. Here, 
we try to locate the threshold value of q for complete 
vacuum instability. 

We perform the defect-type dynamic simulations for 
several q values (10~^ ~ 1.0) to locate Pc- In Fig. 21 
we present simulation results for q — 0.1 at p — 0.999. 
Effective exponents S{t), rj{t) and \/z{t) show upward 
curvatures, which imply that the system is still active 
even at p = 0.999. It suggests that the criticality is 
located at p = 1.000(1) and the vacuum is completely 
unstable. For other nonzero q values, we also find the 
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FIG. 4: Plots of the effective exponents against 1/t at 
p = 0.999 for q — 0.1. Upward curvature of each exponent in- 
dicates that the system is still in the active phase at p = 0.999. 

similar results to those for q — 0.1. We conclude that, 
for any g 7^ 0, the system is always active and only critical 
at Pc = 1.0. As discussed in the previous subsection, this 
result again confirms that the multispecies character is 
relevant (not the symmetry) to vacuum instability. We 
also check the C-dominant regime for g > 1/2 and find a 
similar result. 

To identify the critical behavior, we estimate the ex- 
ponent i>±^ by collapsing the particle density data with 
(3/i'± = 1. Fig.[Slshows the scaling plot for q = 0.1 with 
= 0.5. We estimate v± = 0.50(5) for q = 0.1. For 
other nonzero q values, we also estimate iy± ~ 0.5. 

Although we do not consider the incomplete coupling 
case (cr ^ 0), we expect a nonzero threshold value of a* , 
below which the system is always active. 

IV. CONCLUSION 

In this paper , we study the stability of vacuum in A^ 
coupled DP systems. The interspecies coupling is lin- 
ear, bidirectional, and excitatory. In the coupled contact 
processes, the vacuum is always stable at a sufficiently 
low branching rate for all A^, regardless of the coupling 
strength cr', and the system undergoes DP type absorb- 
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FIG. 5: Data collapse of psL''^^"^ against AL^^"^ with 
j3/v±_ = 1 and v± = 1/2 for g = 0.1. The steady state density 
ps is measured for various system size L (2^ ~ 2^) and A 
(5 X 10"'' ~ 3 X 10"^). 

ing transition into the vacuum. On the other hand, in the 
coupled BAW with one offspring, the vacuum stability is 
quite fragile for > 2 in a strong coupling regime. The 



absence of a spontaneous annihilation process ^ is 
crucial for vacuum instability. 

We find that the vacuum is unstable with respect to 
an arbitrarily small branching rate in a sufficiently strong 
coupling regime (cr/ > cr'*) for > 3. The multispecies 
character is the key element responsible for this vacuum 
instability and the asymmetry in the interspecies cou- 
pling is shown to be irrelevant. The N = 2 case is special. 
The vacuum is always stable like in the coupled contact 
processes, but the HC interaction is relevant to vacuum 
instability in a strong coupling regime. We show that the 
N = 2 BAWi model without HC interaction can reduce 
to the A^ = 1 model even at the full coupling strength, 
which explains the specialty at N = 2. 

Critical behavior near the annihilation fixed point in 
a strong coupling regime can be conjectured by a sim- 
ple argument of balancing two time scales of branching 
and annihilating random walks. Numerical investigations 
confirm our conjecture of /? = 1/2, = 1/2, and = 1, 
which also applies to Ar-BAW2 with static branching and 
HC interaction p^ . 
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